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Abstract 

We study the correspondence between phase-space localization of quan- 
tum (quasi-)energy eigenstates and classical correlation decay, given by 
Ruelle-Pollicott resonances of the Frobenius-Perron operator. It will be 
shown that scarred (quasi-) energy eigenstates are correlated: Pairs of 
eigenstates strongly overlap in phase space (scar in same phase-space re- 
gions) if the difference of their eigenenergies is close to the phase of a 
leading classical resonance. Phase-space localization of quantum states 
will be measured by L 2 norms of their Husimi functions. 



1 Introduction 

Schnirclman's theorem [1] as well as Berry's physical reasoning predict that 
quantum energy eigenfunctions (Wigner or Husimi representation [2, 3]) of sys- 
tems whose classical counterpart is chaotic are uniformly distributed on the 
energy shell. Heller, however, has shown that there exist quantum eigenfunc- 
tions which are strongly localized (scarred) on hyperbolic periodic orbits [4] . At 
present, the general opinion is that scars are exceptional, while the majority con- 
sits of uniformly (phase-space) distributed eigenstates [5]. On the other hand, 
it is well known that there also exist weakly localized eigenfunctions, which 
continuously fill the "gap" between uniformly distributed and strongly localized 
eigenfunctions. Conveniently, instead of single exceptions we here consider lo- 
calization properties of the whole set of eigenfunctions of the system; restricting 
our studies on finite dimensional Hilbert spaces. To that purpose, we introduce 
squared L 2 norms of Husimi eigenfunctions as a measure for phase-space local- 
ization [6]. As it will be shown in the sequel this measure proves amenable to 
semiclassical considerations. 

The classical dynamics of chaotic systems can be described by the time 
evolution of phase-space density functions. The corresponding propagator is 
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the Frobenius- Perron operator [7, 8, 9], where the poles of its resolvent are 
called Ruelle-Pollicott resonances which coincide with decay rates of classical 
correlation functions [10, 11, 12, 13]. The quantum-classical correspondence, in 
particular the influence of Ruelle-Pollicott resonances on the quantum energy 
spectrum, is still a point of interest [14, 15, 16, 17]. We show that phase- 
space overlaps of energy eigenstates turn out Lorentz distributed with respect 
to the differences of their eigenenergies. The Lorentzians are determined by 
Ruelle-Pollicott resonances. In other words, the probability that two eigenstates 
strongly overlap becomes large if the difference of their eigenenergies coincides 
with the position where the Lorentzian is peaked. On the other hand, if a pair of 
eigenstates strongly overlaps (much more than random-matrix theory predicts) 
each of them must be localized, i.e. scarred, in the same phase-space regions. 

We here consider systems with a compact two-dimensional phase space, in 
particular the unit sphere, whereby the Hilbert-space dimension of the quan- 
tum counterpart becomes finite. Periodically driving destroyes integrability in 
general. Moreover, a stroboscopic description leads to a Hamilton map or a 
Floquet operator in the classical or quantum case, respectively. A well known 
representative of such a dynamics is the kicked top [18, 19, 20]. 

2 Kicked Top 

The dynamics of the kicked top is described by a stroboscopic map of an angular 
momentum vector whose length is conserved, For such a dynamics the phase 
spaces is the unit sphere. The classical time evolution is usually described in the 
"Hamilton picture" , whereby the stroboscopic consideration leads to a Hamilton 
map which describes a trajectory after each period, 

(q', P ') = M(q,p). (1) 

Here the primes denote the final position q and momentum p coordinates. On 
the sphere the canonical phase-space coordinates are given by the azimuthal 
and polar angle as q — (p and p = cos 9. 

The quantum time evolution in the Schrodinger picture is generated by a 
Floquet operator F which is built by the components of an angular momentum 
operator J; we choose it as a sequence of rotations about the y and z-&xis 
followed by a nonlinear torsion about the z-axis, 

F = T z (T)R z (a)R y (/3) , 
T z (t) = e-^ J * , 
R z (a) = e~ [aJ * , 

Ry(J3) = e"^ , (2) 

where r is called the torsion strength and a and (3 are rotation angles. The 
dimension of the quantum Hilbert space is TV = 2j + 1 , where j is the quantum 
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angular momentum formally replacing the inverse of Planck's constant, %~ 
(h = 1 in this article). Since F is unitary, it has N orthogonal eigenstates 
with unimodular eigenvalues characterized by eigenphases (quasi-eigenenergies) 
as F n \4>i) = exp(-in0i)|<fc). 

The classical counterpart has in general a mixed phase space. By the choice 
of the parameters as a = f3 = 1 and t — 10 the dynamics becomes strongly 
chaotic. For j = 200 which we use for numerical results stable island are not 
resolved by the Planck cell of size whereby the dynamics looks, from a 
quantum point of view, effectively hyperbolic. 

In order to compare the results of the kicked top with those of random-matrix 
theory (RMT) we here discuss the symmetries of the system. The dynamics 
proves invariant under nonconventional time reversal. In terms of random- 
matrix theory the Floquet matrix belongs to the circular orthogonal ensemble 
(COE), where the coefficients of the eigenvectors can be chosen real in a suitable 
basis [21, 22]. We here expand the Floquet operator in the basis of eigenstates 
of the .z-component of the angular momentum operator, J z \jm) — m\jm). By 
a unitary transformation given by a simple rotation, F' = RFW , the Floquet 
matrix becomes symmetric, F' T — F' , where the eigenvectors become real. 
Here T denotes transposition. This is an important property, since the L 2 
norm of a Husimi function which we will use as a measure for phase-space 
localization is invariant under rotations. Therefore, eigenvectors of the kicked 
top must be compared to real random vectors. The rotation is of form R = 
Rz(%-)Ry{—f)Rz{—f), whereby the transformed Floquet matrix becomes 

F> = R z (§ )R y (-^)R z (-l)T z (r)R z (a)R y (P)R z (l)R y ^ )R z (-§ ) . (3) 

After commutation of the rotation R z (— f ) with the torsion the product 

R z {-*-)R y {P)R z C-) = R y ^)R z (P)R y (-l) (4) 

is a rotation about the rr-axis. Using the latter relation the transformed Floquet 
matrix finally becomes 

F> = (R y {l)R z {t)) T T z {T)R z {a)R y {^)R z {t) , (5) 

which is obviously symmetric. 



3 Frobenius-Perron Operator 

Another way to describe classical time evolution is the "Liouville picture" , as 
the propagation of density in phase space. The corresponding propagator, the 
Frobenius-Perron operator V, is defined through the Hamilton map as 

f n (q,p) = Pf (q,p) = fo(M- n (q,p)) = J dq' dp' f Q (q' , P ')5 2 ((q' ,p')-M- n (q,p)) , 

(6) 



3 



where f(q,p) denotes an arbitrary phase-space density function. Note that the 
Hamilton map is invertible and area preserving. An expectation value of a 
classical observable is given by the phase-space integral 



(A)= dqdpA(q,p)f(q,p) = (A\f) 



(7) 



where we have introduced the Dirac notation. To avoid confusion with quantum 
wave functions we here use round brackets. Note that this notation is generally 
to read as a linear functional, where the density function / belongs to the Banach 
space L 1 and the observable A to the dual space L°°. Furthermore, we suppose 
that both functions are real, otherwise A is to complex conjugate in the integral 
notation (7). 

For classically chaotic systems (we assume purely hyperbolic dynamics) cor- 
relations of observables decay exponentially in time. Due to ergodicity the time 
correlation can be written by a phase-space integral, 



where the time dependence of the observables must be read as A(n) = A(q(n) , p(n)) . 
Here pi denotes the stationary (invariant) density with Vpi = pi, i.e. the con- 
stant on the sphere. The associated stationary eigenvalue is 1 and ensures that 
no probability gets lost, i.e. it preserves the L 1 norm of a density function. We 
may replace the observable B by an initial density function / and further we 
assume that (Alpi) = 0, then the correlation function can be written in terms 
of the Frobenius-Pcrron operator. Finally, we introduce the Ruelle-Pollicott 
resonances A„ which are to identify as decay rates, 



where the a v denote the coefficients of the resonance expansion. We here as- 
sume the simplest case that resonances appear with multiplicity 1, otherwise the 
spectral decomposition of the Frobenius-Perron operator would be given by a so- 
called Jordan block structure, whereby the expansion on the rhs would become 
more complicate. It should be remarked that the decay rate is precisely given 
by the logarithm of \ v which is more convenient to consider in continuous-time 
dynamics. While the A„ are located inside the unit circle on the complex plane, 
the logarithm of \ v are customarily chosen to be in the lower half plane. Due to 
the fact that the Frobenius-Perron operator preserves the positivity of density 
functions the resonances are real or appear as complex pairs. The trace of the 
Frobenius-Perron operator defined through trP™ = J dqdpS 2 ((q,p)~M n (q,p)), 
i.e. setting image and original points in (6) equal, becomes 



C AB (n) = (A(n)B(0)\ Pi ) - (A\ Pi )(B\ Pi ) 



(8) 




(9) 




(10) 
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in terms of the resonances. We here have separated the stationary eigenvalue 1 
from the summation of the resonances. We should carefully distinguish between 
forward and backward time evolution, since we don't expect an increase of 
correlations for the backward time propagation. It has been shown that the 
backward time Frobenius-Perron operator has the same resonances. 

Ruellc-Pollicott resonances arc defined as the poles of the resolvent of V, 

R W = J^p- (11) 

The corresponding "eigenfunctions" are not square-integrable functions like 
those of the quantum propagator, but distributions. It is known that unstable 
manifolds of periodic orbits function as supports of these singular eigenfunc- 
tions. For the backward time evolution stable and unstable manifolds exchange 
their role. 

4 Approximate Resonances 

Supported by arguments of pertubation theory, Weber et al. [23, 24] have shown 
that classical Ruelle-Pollicott resonances of the Frobenius-Perron operator can 
be found by investigating the propagator restricted on different phase-space 
resolutions. Moreover, one finds approximate eigenfunctions which scar along 
unstable manifolds. This important result will be useful for the comparison of 
classical and quantum eigenfunctions. 

We restrict our considerations on the Hilbert space of square integrable func- 
tions L 2 . For the system considered the phase space is the unit sphere, where 
it becomes convenient to use the basis of spherical harmonics. The Frobenius- 
Perron operator becomes an infinite unitary matrix whose unimodular spectrum 
can be separated into a discrete part for integrable components and into a con- 
tinuous part for hyperbolic components of the dynamics. Truncation to an 
M x M matrix corresponding to a restricted phase-space resolution destroys 
unitarity. The spectrum becomes discrete and the eigenvalues are inside or on 
the unit circle. As M increases some eigenvalues prove M-independcnt; these 
are said to be stabilized. Stabilized eigenvalues reflect spectral properties of the 
Frobenius-Perron operator. They arc (almost) unimodular for integrable com- 
ponents or stable islands. In contrast, eigenvalues which are stabilized inside the 
unit circle reflect Ruelle-Pollicott resonances. Non-stabilized eigenvalues have 
typically smaller moduli than the stabilized ones. They change their positions 
as M increases till they reach positions of resonances where they can settle for 
good. We also know that the eigenfunctions corresponding to stabilized eigen- 
values are localized either on tori for unimodular eigenvalues or around unstable 
manifolds for non- unimodular eigenvalues. We expect that the latter eigenfunc- 
tions converge weakly to singular resonance "eigenfunctions" . We will call them 
approximate resonance eigenfunctions in the sequel. 
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Fig. 1 (a) shows the eigenvalues of the truncated Frobcnius-Perron matrix 
in the complex plane. Here the dimension is M = (7 max + l) 2 , where Z max is 
the maximal total angular momentum of the spherical harmonics in which the 
density functions are expanded. In (b) we see a grey-scale shaded plot of the 
eigenvalue density calculated from truncated matrices (7 max = 20, 21, ... , 70). 
Dark spots corresponding to large amplitudes of the density indicate resonance 
positions, since resolution independent eigenvalues highly increase the density 
through accumulation. A comparison of both shows that some eigenvalues in 
(a) reflect resonance positions (see also Tab. 1). In Fig. 2 (a) the modulus of 
the approximate eigenfunction is plotted in phase space, where the dark regions 
belong to large moduli of the complex-valued function. The comparison with 
the unstable manifolds (b) of a weakly unstable period-4 orbit shows that the 
approximate eigenfunction scars along these manifolds. In (c) and (d) we see 
the scarring of the eigenfunction of backward time propagation (same resonance 
as in (a)) along the stable manifolds. 



5 Coherent-State Representation 

To represent quantum operators in a way close the corresponding classical ob- 
servables it is convenient to start from coherent states. For the SU(2) group 
coherent states can be generated by a rotation of the state \j,m — j) as 
\0,<p) = R(9, ip)\jj}- In the \jm) basis these coherent states are given by 

M=(l+(tan§) 2 ) 3 £ JUQ (tan§ e ^y'" ro |im>. (12) 

\ / 771= — j V ^' ' 

Since the set of the coherent states is overcomplete, there are several ways 
describing quantum operators with coherent states [25, 26, 27]. On the one 
hand, one can use the so-called P function, defined as the weight of coherent- 
state projectors in the continuous mixture 

A=£-JtilP A (0,<p)\6,<p){6,<p\. (13) 

The integral is over the unit sphere, where dQ = d(pd8 sin#. On the other hand, 
we have the so-called Q function, the coherent state expectation value of an 
operator, Qa(9, <p) = {9, y>\A\6, ip). It is important that Q is a smooth function 
on phase space, while P can strongly oscillate, particularly in the shortest wave 
lengths. However, in contrast to the coherent states of Weyl groups, P functions 
always exist. Both functions can be expanded in terms of spherical harmonics 
[28], 

2j I 

Q A {9, V ) -EE lUA)Yr(6,<p) 

1=0 m=-l 
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2j I 

P A (6,<p) = Yl E Pim(A)Yr(6,v), (14) 
;=o m=-l 

where the summations break off at the total angular momentum I — 2j. The 
Hilbert space of phase-space functions must not be confused with the Hilbert 
space of quantum wave functions; we therefore use round brackets for the scalar 
product already introduced in (7), (f\g) = J dSl /* (8, (fi)g(9, <p). Although the 
P functions tend to oscillate more strongly than the Q functions, the expansion 
coefficients of the P and Q functions corresponding to the same operator con- 
verge to one another in the classical limit, pi m (A) r ^^q lm (A) for fixed l,m. It 
is easy to see that the trace of an operator product can be written as the scalar 
product of P and Q as tiA^B — j^(Pa\Qb)- In particular, if the set {|</>fc)} 
of wave functions form an orthogonal basis of the quantum Hilbert space, then 
the P and Q functions of ket-bras Pik = P\4>i)(<j> k \ and Qik = Q\<t>i)(4> k \ generate 
biorthonormal sets in the Hilbert space of phase-space functions 

N 

-r(Pik\Qt>k') = tr(|0feX<fc||&'X0k'l) = hk'Sw • (15) 

The Q function of a density operator Q p is also called Husimi function. If 
the density operator is a projector of form \ip)(il>\ the corresponding Husimi 
function is a phase-space representation of a quantum wave function. We call 
the Qik = Q\<t>i){4> k \ Husimi eigenfunction if the \<pi) denote Floquet eigenstates. 
The latter notation becomes obvious in the next section. We further distinguish 
between diagonal Husimi eigenfunctions Qkk and skew ones Qik with i ^ k. 



6 Husimi Propagator 

The Husimi propagator T is defined through the time evolution of a quantum 
density operator p as 

Q p{ n)(0,y) - (6,p\F n p(0)(Fi) n \0,v) =F n Q p(0) {e,<p). (16) 

Using the Floquet eigenstates, F\4>i) = e -1 ^!^), the Husimi eigenfunctions 
are easily calculated as Qik = (0 ', (p\<j>i) (4>k\& ', <p) ■ The Husimi propagator thus 
has iV 2 unimodular eigenvalues whose phases are differences of the Floquet 
eigenphases. There is an AT- fold degeneracy of the eigenvalue 1 corresponding 
to the diagonal Husimi eigenfunctions Qkk which are real and normalized as 
J d^lQkk — Air/N. All other (skew) Husimi eigenfunctions with i ^ k are 
complex and their phase space-integral vanishes. In basis of spherical harmonics 
the Husimi propagator becomes an N 2 x TV 2 matrix (Husimi matrix). The 
diagonal representation of the Husimi propagator is simply given as 

F n - E lfeK in( ^~ 0fc) ( p ik\ ■ (17) 

i=0 k=0 
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The Husimi spectral density is identified as the density-density correlation func- 
tion with respect to the Floquet eigenphases, 

/>27T 

C(u>)= d#(" + = -(&-&)) , (18) 

where p(<fi) = J2 k 8{<j) — <j)k)- Some authors prefer the normalized density p(w) = 



7 L 2 Norms of Husimi Eigenfunctions 

The squared L 2 norm of a diagonal Husimi eigenfunction, 

||Q fcfc || 2 = J dSl\(0,tp\<t>k)\* , (19) 

is the inverse participation ratio (IPR) with respect to coherent states (phase- 
space distribution) . It becomes large if the Husimi function is strongly localized 
in phase space, say scarred on periodic orbits. On the other hand, the squared 
L 2 norm of a skew Husimi eigenfunction can be understood as the overlap of 
two diagonal Husimi functions on phase space, 

dn\Q lk \ 2 = J dil\{9,<p\<l>i)\ 2 \{0M<i>k)\ 2 ■ (20) 
From the Schwarz inequality, 

||Oifc|| 2 <||Q«||||Qfcfc||, (21) 

it becomes obvious that for large values of ||<2i/c|| 2 both diagonal Husimi eigen- 
functions must be localized in the same phase-space regions. 

We illustrate two examples. A constant function on the sphere is of course a 
uniformly distributed function. But note that it is not a Husimi eigenfunction, 
since the corresponding density operator is of form p = j^l. Using the normal- 
ization £j JdflQp = 1 the Husimi function becomes For the squared L 2 
norm we find 

4-TT 

WQp\\ 2 = (22) 

for uniform function. In comparison, the random-matrix averaged squared L 2 
norms of diagonal Husimi eigenfunctions is larger by a factor 2 (see Sec. 8) which 
can be explained by quantum fluctuations. A most strongly localized Husimi 
function, in contrast, corresponds to a density oparator which is a coherent-state 
projector. Due to the invariance under rotations the L 2 norm is the same for 
all coherent-state Husimi functions. Using the coherent-state projector \jj)(jj\ 
and (27) one easily finds for the squared L? norm 

2-7T 

IIOi«x«ill 2 = (23) 
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for strongly localized Husimi function. It becomes obvious that squared L 2 
norms of most strongly localized and uniformly distributed Husimi functions 
differ by a factor of order N. 

We here show the calculation of the L 2 norm from vector coefficients in the 
\jm) basis (c^ = (jm\4>k))- Introducing the decomposition of unity in terms of 
the angular momentum states, 1 = X^m=-j b m )O m lj the L 2 norm of a skew 
Husimi eigenfunction becomes a four-fold sum, 

/d^,^><<fel^>l 2 = £ C(4J*4,(0* 

mi ...7714 

x J dfl{9,ip\jm 1 ){jm 2 \e,ip){e,ip\jm 3 ){jm i \9,ip) . (24) 



From the definition (12) of the coherent states we find, using the relations 
tan % = ±=2** = jj^-z for < 9< tt, 

2 sin $ 1+cos 



(jra|<M = 2~ J 'y ^. J* m )(! - 008 0)^(1 +cos0) i ^e i ( J '- m ) ¥ ' . (25) 

After a few further steps the phase-space integral f dQ = f£ d6 sin 6* J" Q w dtp 
leads to a Kronecker 5 and to a beta function [29] for the (p- and ^-dependent 
part, respectively. We thus obtain for the coefficients of the summation (24) 

J dfi f\jmi){jm 2 \0, ip) (9, (p\jm 3 ) (jm 4 \9, if) 

r( ^z^ + i)r(^±& + i) 



= 47T 



r(4j + 2) 



x<J(m 1 -m 2 + m 3 -m4)W( . ^ )( . 2j )( . 2j )( . 2 j ] (26) 



v .j - mi J \j - m 2 J \j - m 3 y \j - m 4/ 
We make use of the argument of the Kronecker S and get for the L 2 norm 
f Ao\ia \m\ia. \a \i2 a (2j-mi -m 3 )!(2j + mi +m 3 )! 

y dft|<0,¥>|<fc)<<fo|0,p)| = 4tt 2^ (iJTT)! 



2j W 2j \ / 2j \ / 2j 
j - mj \j - m 2 ) \j - m 3 J \j - m 4 

x 5{mi-m 2 +m 3 -m i )c\ ni {c l m2 )*c^ 3 {c^ li )* . 

(27) 

The latter formular will be used for the calculation of Husimi L 2 norms of 
random vectors in the next section. 
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8 Random-Matrix Average of L 2 Norms 



In order to compare the results of the kicked top with those of random-matrix 
theory we now consider diagonal Husimi eigenfunctions and replace the coef- 
ficients by real or complex random numbers due to the orthogonal (COE) or 
unitary ensemble (CUE), respectively [30]. It is easy to see from (19,20) that the 
L 2 norm of a Husimi function is invariant under rotations. As it has been shown 
in Sec. 2 the eigenvectors of the kicked top must be compared to real random 
vectors. The only correlation between the coefficients is the normalization of 
the vector, ^2 m \c m \ 2 = 1 (we drop the upper index in the following). Therefore 
we neglegt all terms containing random phases and keep the contributing terms 
with mi = 77i2 A ?TJ3 = TO4 or mi = m4 A m,2 = 777,3. By the choice of real coeffi- 
cients one has the further possibility mi = m3 A m-2 = 777.4. We may abbreviate 
the coefficients of the summation (26) as /(mi, . . . ,m 4 ) for a moment to find 
out the contributing terms, 

(IIQfefell 2 ) = f( m i,---, m 4) S { m i- m 2+m 3 -Tn 4 )(c mi c* m2 c m3 c* mi ) 

mi.. .7714 

= f( m ii • ■ ■ > m 4) S(m 1 -m 2 -\-m 3 -m4) 

mi.. .7714 

X ^^77117712^77137714(1 ^771 17713 ) H~ ^77117714 ^7712777,3 (l ^777-17773 ) 

H" {^77117773 ^77727714 (1 _ ^^2)} ( | C m | | C 7l | ) + ^m;m fc ( I C 77l | ) j ■ 



(28) 

The diagonal part S^ n . mk = S mim2 S m2m3 S m3m4 is separated, since there appears 
the average of |c m | 4 . Due to the symmetry of /(mi, . . . , 777,4) the first two terms 
in the bracket give the same contribution; to that end, one executes the sum- 
mations over m 2 and m 4 . The contribution coming from real coefficients {. . .} 
vanishes; here one summates over 7713 and 7714. Finally, one finds for the averaged 
L 2 norm 



(WQkkW 2 ) 



47T 



4j + 1 

3 

+ E 



2 JC 

2j - 2m 



4.7 

2j — m — n 



2j 

3 - m 



2j \f 2j 
j-mj \j - n 



<|c™| 4 ) 



Next we need the averages of products, 



CUE 



(| c m| \ c n\ ) \ [ 1 , 1. ,r 

I N(N+2) ' ^ W£j 



I l 2 l l 2 \ 
I Cm I I Cn I / 

(29) 
(30) 
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- , CUE 

N(N+2) 



<w 4 > = \ ^ ;coe ■ ( 31 ) 



The averages are easily calculated from the probability distribution given in [30] . 
It should be remarked here that semiclassical corrections are of next to leading 
order of TV -1 . Therefore we do not neglect this order in our RMT results. 
Let us consider first the CUE. Since (|c m | 4 ) = 2(| 

Cm | | C n | 2 ), we can com- 
plete the summations over the diagonal (m = n) and off-diagonal (m ^ n) part. 
Thus we have for the CUE 

(HOfefeii 2 ) = N{N+i)4j + i J2_ _ ( 2j _ m _ n ) (j _ m ) (j _ n ) • 

(32) 

It might not be easy to see that this results to (||Qfcfc|| 2 ) = jv(jv+i) ■ To that 
end we present an alternative way to calculate the averaged L 2 norm which is 
unfortunately not applicable for the COE [6]. An arbitrary complex random 
vector can be written as = U\jj), where U is a unitary random matrix. A 
coherent state is generated through a rotation operator as \0,ip) — R(0,<p)\jj). 
We thus have (0,<p\ip) = (jj\R^U\jj) . The product R)U = U defines a new 
unitary random matrix, where the Haar measure keeps unchanged d[i(U) = 
d(j,(U). Therefore the averaged L 2 norm becomes 

(WQkkW 2 ) = J d»(u) J dn\(jj\rfu\jj)\ 4 
dn [ d(,(u)\(jj\u\jj)\ 4 



= 4^(|c m | 4 ). (33) 

Note that for the COE, in contrast, one deals with real vectors which generally 
become complex after rotating. 

For the COE we again complete the diagonal and off-diagonal summations 
in (29), whereby a further contribution remains, because of the factor 3 of the 
fourth moment (31), 

4tt ^ ( 4j \"7 2j y 1 (4 J (2j)!) 2 4tt 

7 ^ 



4j + l ^_\2j-2m) \j-mj N(N + 2) (4j + l)(4j)! N(N + 2) 

m — j 



The latter equation can be calculated as follows [32]. We rewrite 



(34) 



V { A] W 2j V (2j)! 2 yV2fc\/4j-2fc\ 

2. { 2 j - 2m) {j -m) (4j)! 2^\k)\2j-k) 1 > 
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with k = m + j. We now show that 



r 2k\ (2n - 2k 



k=0 



S(7 „-D =4 "- < 36 » 



The generating function of ( 2 fc fc ) is 



which is easily seen from Taylor expansion, 



1 -tvv 



k\dx k[ ' 



fc!2 1+ 2 - fc -2 4 (38) 



z=0 

Squaring the generating function gives 

^=-f:f;(t)(?)--£-±(?)C::r). « 

fe=o ;=o v / \ / „ = o fe = o v / \ / 

where n = k + I. Comparing the latter equation with the geometric series 

_ = £y4» (40) 

finishes the proof. 

The aforementioned contribution (34) becomes more familiar by the approx- 
imation _ 

(4^(2j)!) 2 _V^r(2. ? + l) _1 (41) 



(4j + l)(4j)! 2 T(2j + 3/2) 2^2j + \- 
The averaged L 2 norms finally become 

<««-^=^{ A(^vm ^ • (42) 

As an interesting point to remark here the COE cigenstates are somewhat more 
localized than the CUE eigenstates, since their averaged L 2 norm contains a 
further contribution of order . But in contrast to the fourth moment (31) - 
it might be understood as the averaged IPR with respect to the \jm) basis - the 
difference vanishes in the classical limit N — > oo. However, for both ensembles 
the averaged squared L? norms are roughly twice larger than for a constant 
distribution on the sphere. 
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For the skew Husimi eigenfunctions the averaged L 2 norms can be calculated 
from the relation 



4tt = J dci\(e, v \e,^)\ 2 = J dn l^ne^i^iA = EH^H 2 +Ell^" 2 



which leads to 



ik | 

i^k 

(43) 

(WO H2X ^-N(\\Qkk\?) 
Introducing the averages of the diagonal ones we find for CUE 

For the COE one obtaines the same result up to an order of N~i, 

/,,^ „2\ 47r ^ (, 1 47r N «/»r ^ 

(46) 



9 Comparison of Quantum and Classical Eigen- 
functions 

It has been shown in [31] that quantum quasiprobability propagation looks 
classical if phase-space resolution is blurred such that Planck cells are far from 
being resolved. As a result, if we truncate the Husimi matrix to coarse phase- 
space resolution, then it becomes almost equal to the truncated Frobenius- 
Perron operator, TTT N —^TVT, where 

'max I 

Z=0 m=—l 

denotes the truncation projector restricting the "classical" Hilbert space dimen- 
sion to M = (l max + l) 2 . In the following we choose M -C iV 2 , and thus have 
TTT k> TTT. Again we consider the classical propagator. Let \w) — T\w) be 
an approximate resonance cigenfunction with stabilized eigenvalue A of TVT . 
Since the stabilized eigenvalues reflect spectral properties of the non-truncated 
Frobenius-Perron operator, A" must be a stabilized eigenvalue oiTT jn T, at least 
if M is chosen large enough. This property does not hold for non-stabilized 
eigenvalues, since TV a T ^ (TTT) 12 . Choosing the approximate resonance 
eigenfunction L 2 normalized, the return probability becomes (u>|"P|u;) = A" 
for small n (see Fig. 3). Note that the latter property makes sense only if we 
consider approximate resonance eigenfunctions, because the singular resonance 
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eigenfunctions are not square integrablc. Fig. 3 shows the return probability 
of the approximate eigenfunction from stabilized eigenvalue X v w — i0.75 (No 
4 in Tab. 1). As well the moduli as the phases coincide with A" for about 20 
iterations. Beyond this time quantum fluctuations become visible. 

Since TTT w TVT, we expect TT n T w TV n T for small n. Now we can 
replace the Husimi propagator by its diagonal representation (17), 

A™ « (u^H = ^ E (MQik)e- in{ ^- M (Piklw) . (48) 

The return probability becomes a double sum of overlaps of quantum and clas- 
sical eigenfunctions. For large N the Husimi cigcnphases are quite dense in 
the interval [0, 2ir). Outside an interval around the Husimi eigenphase to = 
wherein level repulsion of Floquet eigenphases become perceptible the spectral 
density of differences — <f)k is almost constant, N 2 /2ir. It is convenient to re- 
place the sum by a continuous integral over the Husimi eigenphases, where the 
overlaps can be replaced by a continuous function as a smoothed distribution 
of overlaps, 

A " ~ dw HQ*R AW (w)e- inu ■ (49) 

o 

In the classical limit we let first N — > oo and then M — > oo. In this limit 
the foregoing result becomes valid for all n. The smoothed overlaps are given 
by the inverse Fourier transform of A" as 



{w\Qik){Pik\w) M^]v i + |^-2|A|c s(c-argA) ' (50) 

which is a periodic "Lorentzian" displaced by the phase of A and of width — In | A| 
[33] (see Fig. 4). 

The smoothing is done by a convolution with a sine function which naturally 
results from a truncated Fourier transform. Here the time restriction is J\f < 
n < J\f, where TV w 20 corresponding to the validity of the return probability 
of the approximate eigenfunction (Fig. 4). For simplification the integral is 
restricted on the interval between the first zeros of the sine function, 

/ (w) cx / dui — f(ui -uj). (51) 

7T 

A simple argument connects quantum scars with Rucllc-Pollicott resonances. 
Classical resonance eigenfunctions are scarred along unstable manifolds (or sta- 
ble manifolds for backward time propagation). Quantum eigenfunctions which 
strongly overlap with resonance eigenfunctions have to be scarred as well. This is 
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an important result, because it explains scarring of quantum eigenfunctions not 
only on periodic orbits, but also along stable and unstable manifolds. This has 
been observed for the kicked top [34] . For instance, we consider the skew Husimi 
eigenfunction which shows the largest overlap with the classical approximate 
eigenfunction considered before. The corresponding diagonal Husimi eigenfunc- 
tions (No 11 and No 30 in Tab. 3) plotted in Fig. 5 are scarred in the same phase- 
space regions, where the the approximate resonance eigenfunctions are scarred 
(Fig. 2), while the difference of their Floquet eigenphases (</> 30 — <pu = —1.586) 
is close to the phase of the resonance (arg A4 = — §)■ 



10 Resonance Corrections of Averaged Phase- 
Space Overlaps 

In the foregoing section a qualitative explenation of the connection between 
resonances and quantum eigenfunctions was given, but now we are interested in 
more quantitative results. To that end, we consider transition rates of coherent 
states in the classical limit 

^\{0,v\F n \6',</)\\-^ i 5((q,p)-M n tf,i/)) . (52) 

The latter relation becomes obvious if one suggests that coherent states are wave 
functions most strongly localized on phase-space points. We now consider the 
return probability and integrate over the phase space, 



N 

47T 



<m\(6,<p\F n \e,<p)\ 2 N -^> dQ8((q,p)-M n (q,p)) , (53) 



where the rhs is the trace of the Frobenius-Perron operator [37]. We remark 
here that the integral on the rhs leads to a sum of contributions from periodic 
orbits, which is an important connection between scars on periodic orbits and 
the results of our paper. On the one hand, periodic orbits, in particular the 
weakly unstable ones, contribute to the trace of the Frobenius-Perron operator, 
i.e. influence the resonances. On the other hand, scars typically appear around 
weakly unstable periodic orbits. On the lhs of (53) we introduce the diagonal 
representation of the Floquet operator and get 



52\{0M<t>k)\ 2 e 



f J2WQik\\ 2 e- in ^-^. (54) 



47T 

ik 
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Fourier transformation of the latter expression leads to a sum of S functions 
weighted by L 2 norms, 



E ^Eii^ii 2e " inWi "** ) = Eii^*ii 25 ( w -^-^))- ( 55 ) 

n— — oo ik ik 

Due to the arguments of the foregoing section we expect that for finite N the 
relation (52) is valid for finite times \n\ < N. Validity of semiclassical methods 
is guarantied for times up to the Ehrenfest time, where the number of fixpoints 
coincides with the number of Planck cells. Thus we identify TV as Ehrenfest 
time. The truncated Fourier transform leads to a sum of smoothed 5 functions 
(51). Using (53) we get 

ra^n = % e trvn ^ ( 56 ) 

which we may call smoothed L 2 norms. The next step is to drop the station- 
ary eigenvalue 1 in the traces of the Frobenius-Perron operator. The Fourier 
transform of this eigenvalue leads to a 5 function in the limit J\f — > oo which 
is not a point of interest here. In the Husimi representation we identify the 
eigenvalue 1 as the sum of squared L 1 norms of Husimi eigenfunctions. This 
is easily seen from the Husimi matrix in basis of spherical harmonics. In the 
first row and column is only one non- vanishing matrix clement (} / o°|.7 r |Y ) = 1. 
Note that Y® is a constant function on the sphere and therefore proportional 
to the stationary density, i.e. it is the eigenfunction from eigenvalue 1. Intro- 
ducing the diagonal representation of the Husimi propagator (17) one identifies 
£Ek (y \Qkk)(Pkk\Y Q °) = 1. Note that (Y°\Q lk ) = (P lk \Y°) = ^S ik , where 
(Xo\Qkk) HQfefclli- For n = the trace of the Frobenius-Perron operator is 
not defined. The integral on the lhs of (53), however, is defined and gives the 
leading order contribution N. We replace the traces by sums of the Ruclle- 
Pollicott resonances (10) and make use of the symmetry tr"P~™ = tr"P", 

P^F A >) = + 1 J2 E A " cos ™ ■ ( 5? ) 

n— 1 v 

Note that the eigenvalue 1 is also dropped in the leading order term. 

To get from the smoothed L 2 norms ||Qifc|| 2 (uj) toamean value (||Qifc|| 2 ) Aw ( 
we have to divide it by the smoothed level density of the Husimi spectrum. The 
density of the Husimi spectrum is identified as the density-density correlation 
function with respect to the Floquet eigenphases, 



p H {uj) = J2 S ( UJ - {4>i - 4>k)) ■ (58) 
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This spectral density can be calculated from the time Fourier transformation of 
the form factor. The smoothed density is given by a truncated Fourier trans- 
formation as 

N 2 / 2 M \ 

where we have separated the leading order term (n = 0). It is known that the 
form factor is small as the time n is small ((|trF™| 2 ) = n or w 2n for CUE or 
COE, respectively). Thus, the summation up to the Ehrenfest time which is 
much smaller than N becomes negligible and the smoothed spectral density is 
nearly constant. Before we write down the final result, we consider the truncated 
Fourier transformation of the resonances. Since the moduli of the resonances 
are smaller 1, the summations in (57) converge quickly such that we can replace 
the Ehrenfest time Af by oo, 

<iifeii 2 )M = ^(^ + ^EE A > os -) • (so) 

\ V 71=1 / 

The constant term coincides with the RMT result (46), where = + 
0(N~ 2 ). The resonances lead to an (alias Ti) correction in form of overlapping 
Lorentz distributions (see (50)). In the classical limit the resonance corrections 
vanish as N goes to infinity. However, for finite dimension we see non-universal 
corrections which are related to the chaoticity of the system. If, for instance, 
the classical dynamics is strongly chaotic such that all correlations vanish after 
one iteration, the resonances are close to the origin and the averaged L 2 norms 
show no deviations from RMT result. 

Before we come to numerical results we should discuss some prelimaries. The 
kicked top is known to have a mixed phase space. Although elliptic islands of 
stable periodic orbits are much smaller than the Planck cell, bifurcations can be 
responsible for further localization phenomena which are sometimes called super 
scars [35]. But in contrast to the resonances, bifurcations strongly influence 
spectral correlations. That means that peaks resulting from bifurcating orbits 
are higher for the smoothed L 2 norms than for the averaged L 2 norms. Thus, 
we are able to distinguish between localization phenomena of resonances and 
bifurcations. The smoothing is again done by a convolution with a sine function 
(51). From the contributions of the diagonal Husimi eigenfunctions we have 
neglegted the squared L 1 norms which correspond to the stationary eigenvalue. 

Fig. 6 shows (a) the smoothed L 2 norms, (b) the mean L 2 norms, and (c) the 
smoothed spectral density. In (a) we see a couple of peaks at the Husimi eigen- 
phascs u = 0, 7r, ±-|, and i^p- In (c) are remarkable peaks at the positions 
u) = ± 2 j- and after division we see in (b) that these peaks are suppressed, while 
the other peaks still have the same magnitude. Comparison with the semiclas- 
sical prediction is shown in Fig. 7. Due to the fact that stabilized eigenvalues 
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representing the resonances are of largest moduli we expect that the semiclassical 
prediction is almost independent of the set of eigenvalues as well as all stabilized 
ones are taken into account. In (a) the semiclassical prediction is computed from 
(i) eight stabilized eigenvalues, (ii) eigenvalues of modulus larger 0.45, and (iii) 
almost all eigenvalues of the truncated Frobenius-Perron matrix. Comparison 
of (i)-(iii) shows that the semiclassical prediction is mainly influenced from a 
few classical eigenvalues of large moduli which represent the resonances. In (b) 
we compare the semiclassical prediction (iii) with the quantum result. In par- 
ticular for the peaks we find a very good agreement. This result shows that 
the probability to find strongly overlapping eigenfunctions becomes large if the 
differences of their eigenphases coincides with the phase of a leading resonance, 
i.e. resonance of large modulus. In comparison with the background the peaks 
are small (a few percent). However, we show in the next section that strongly 
scarred eigenfunctions are mainly responsible for the peaks. The eigenvalues of 
the Frobenius-Perron matrix (c) are plotted logarithmically (In |A| versus arg A). 
These are easily associated with the peaks in (a). It should be remarked that 
there is no eigenvalue of large modulus which corresponds to the small peaks at 
u = ±^f in Fig. 6 (a). 

11 Scarred Eigenstates of the Kicked Top 

In this section we consider single eigenfunctions and verify the statement that 
phase differences of strongly overlapping eigenfunctions coincide with phases 
of leading resonances. Due to the Schwarz inequality we first check that skew 
Husimi eigenfunctions with large L 2 norms are composed by scarred diagonal 
Husimi eigenfunctions. In Tab. 2 we find 25 L 2 norms of most strongly lo- 
calized skew Husimi eigenfunctions, their eigenphases, and the corresponding 
diagonal Husimi eigenfunctions (the numbers correspond to the enumeration in 
Tab. 3). Due to the symmetry Qn- = Q* ki we have restricted the eigenphases 
as < u < ir. The 32 most strongly localized diagonal Husimi eigenfunctions 
are presented in Tab. 3. Comparison of both tables proves that all skew Husimi 
eigenfunctions considered are composed from at least one scarred diagonal eigen- 
function. In Fig. 8 (a) L 2 norms of all diagonal Husimi eigenfunctions are shown, 
while in (b) we see 530 (of total 1.6 x 10 5 ) L 2 norms of most strongly localized 
skew Husimi eigenfunctions. Due to our foregoing results localized skew Husimi 
eigenfunctions appear frequently around the resonance phases. Interestingly, 
there are only two remarkable eigenfunctions (No 11 and 12 in Tab. 2) around 
the Husimi eigenphases lo = ±^p- Moreover, in both cases one of the under- 
lying diagonal Husimi eigenfunctions is No 16 in Tab. 3. Further investigation 
of this eigenfunction has shown that it is strongly scarred on two bifurcating 
orbits of periods 1 and 3. It seems that we found a super scar corresponding to 
a period-tripling bifurcation. 
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12 Conclusion 



In conclusion, phase-space localization of quantum (quasi-)eigenenergy func- 
tions, say scarring, is not only explained by periodic orbits, but also by Ruelle- 
Pollicott resonances and their corresponding resonance eigenfunctions. In par- 
ticular, we found the interesting result that quantum Floquct eigenfunctions 
are pairwise localized in the same phase-space regions if the difference of their 
(quasi-)eigenenergies coincides with the phase of a leading resonance, i.e. reso- 
nance close to the unit circle. But note that this is a statistical statement which 
does not make a prediction for individual eigenstates. Moreover, we can not 
determine if there are either a few strongly scarred or many weakly localized 
eigenfunctions. However, the semiclassical prediction of the averaged L? norms 
is in a good agreement with numerical results. 

The correspondence between scars around periodic orbits described by Heller 
and the results of this paper might be understood as follows: resonances can be 
computed by a so-called cycle expansion, where resonances appear as roots of a 
polynomial whose coefficients are calculated from contributions of short periodic 
orbits (pseudo orbits) [36, 37]. On the one hand, scars typically appear around 
weakly unstable periodic orbits. On the other hand, these weakly unstable 
orbits mainly induce the cycle expansion. 

Although the kicked top has a mixed phase space, localization effects of 
stable orbits or bifurcations can be neglected if such phase-space structures 
are not resolved by the Planck cell. Further investigations are needed for the 
understanding of the so-called super scars which are related to bifurcations. 

For fruitful discussions the author thanks Fritz Haake, Karol Zyczkowski, 
Sven Gnutzmann, Joachim Weber, Petr Braun and Pierre Gaspard. 

We gratefully acknowledge support by the Sonderforschungsbereich Unord- 
nung und grofic Fluktuationen of the Deutsche Forschungsgemeinschaft. 
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No 


RcA 


ImA 


RcAj, 


ImA, y 


1 


0.811497 





0.81 





2 


0.748029 





0.75 





3 


-0.734887 





-0.745 





4 


0.003495 


-0.733015 





-0.75 


5 


0.003495 


0.733015 





0.75 


6 


0.672874 





0.67 





7 


-0.669930 





-0.68 





8 


-0.611201 





-0.63 






Table 1: Eigenvalues of the truncated Frobenius-Perron matrix (left column) 
which represent the expected resonance positions (right column). 



No 


\4>i - <t>k] 


^rllQifcll 2 


No 


No 


1 


0.02293 


4.87904 


9 


8 


2 


0.02995 


2.17252 


4 


2 


3 


0.03690 


2.12795 




8 


4 


0.12504 


2.26266 


4 


1 


5 


0.66008 


2.53247 


8 


4 


6 


0.69002 


2.14943 


8 


2 


7 


0.95140 


2.07455 


32 


8 


8 


1.58566 


3.28086 


30 


11 


9 


1.61324 


2.07295 


30 




10 


1.63071 


2.31385 


30 




11 


1.93314 


2.77164 


16 


9 


12 


1.95607 


3.47930 


16 


8 


13 


2.12224 


2.03752 


16 





No 


\4>i ~ 4>k\ 


^rllQjfcll 2 


No 


No 


14 


2.41853 


2.01226 


21 




15 


2.61615 


2.08038 


16 


4 


16 


2.64610 


2.04506 


16 


2 


17 


2.92433 


2.26189 


21 


5 


18 


2.97684 


2.21163 




8 


19 


2.99690 


2.05725 




9 


20 


3.01983 


2.73944 




8 


21 


3.10287 


2.10537 


23 


9 


22 


3.11230 


3.05349 


24 


8 


23 


3.11857 


2.86021 


17 


30 


24 


3.12580 


2.55565 


23 


8 


25 


3.13523 


2.48017 


24 


9 



Table 2: Husimi eigenphases, L 2 norms of most strongly localized skew Husimi 
eigenfunctions, and the corresponding diagonal Husimi eigenfunctions (enumer- 
ation in Table 3). Due to the rescaling the RMT- average is nearly 1. 
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No 


(t>k 


^rllOfefcll 2 


1 


0.11280 


3.16342 


2 


0.20789 


2.91329 


3 


0.23200 


2.63425 


4 


0.23784 


3.44206 


5 


0.39507 


2.88719 


6 


0.46842 


2.98415 


7 


0.64437 


2.98154 


8 


0.89792 


7.19529 


9 


0.92085 


4.32135 


10 


1.41602 


2.67399 


11 


1.42231 


2.92272 


12 


1.72257 


2.79565 


13 


1.82172 


2.86009 


14 


1.98572 


3.39027 


15 


2.33226 


2.69246 


16 


2.85399 


4.36732 



No 




^rllQfefcll 2 


17 


2.95522 


2.69166 


18 


3.16445 


2.73181 


19 


3.34145 


2.64167 


20 


3.62419 


2.80308 


21 


3.75393 


3.37501 


22 


3.84585 


2.75931 


23 


4.02372 


2.70190 


24 


4.06881 


2.71314 


25 


4.15409 


2.68584 


26 


4.36618 


2.79003 


27 


4.91979 


2.81994 


28 


5.01485 


2.85668 


29 


6.11281 


2.76855 


30 


6.11984 


6.27044 


31 


6.21927 


2.79705 


32 


6.22970 


2.86170 



Table 3: Floquet eigenphases and 1? norms of most strongly localized diagonal 
Husimi eigenfunctions. Note that the RMT-average is about 2.034 
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-1 1 

Fie A 

Figure 1: (a) Eigenvalues of the truncated Frobenius-Perron matrix (7 max = 32). 
(b) Eigenvalue density computed from Frobenius-Perron matrices with Z max = 
20, 21, . . . 70. The centered disc is not shown, because of the increasing density 
at the origin (see (a)). 
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Figure 2: (a) Grey-shade coded phase-space plot of the modulus of the approxi- 
mate resonance eigenfunction (7 max = 32) from A w — i0.75 (No 4 in Tab. 1). (b) 
Unstable manifold of a weakly unstable period-4 orbit (spots) supporting the 
resonance eigenfunction. (c) Approximate resonance eigenfunction of backward 
time evolution from same resonance as in (a), (d) Stable manifold of the same 
orbit as in (b). 
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Figure 3: Return probability of the approximate resonance cigcnfunction from 
A w — i0.75 (No 4 in Tab. 1; see also Fig. 2 (a)): (a) the amplitude shows for 
small n an exponential decay of (0.75)™ (dashed), while the phase (b) evolves 
like — n^. 



i_ 0.00012 
8e-05 








Figure 4: (a) Overlaps of Husimi eigenfunctions with the approximate resonance 
eigenfunction from A w — i0.75 (No 4 in Tab. 1). (b) Smoothed overlaps (solid) 
in comparison with the Lorentzian distribution (dashed) corresponding to A 4 = 
-i0.75 (see text). 
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Figure 5: Grey-shade coded phase-space plots of the diagonal Husimi eigenfunc- 
tions (a) No 30 and (b) No 11 of Tab. 3. The Husimi function in (a) is strongly 
scarred in phase-space regions, where the approximate resonance eigenfunctions 
are scarred (Fig. 2). The Husimi function in (b) is also scarred in the same 
phase-space regions, but shows more structures all over the phase space than 
the function in (a). 
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1 2 3 4 5 6 
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Figure 6: (a) Smoothed and (b) averaged 1? norms 
tions. (c) smoothed (Husimi) spectral density. 
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Figure 7: (a) Semiclassical predictions computed from (i) 8 stabilized eigen- 
values (dotted), (ii) eigenvalues with |A| > 0.45 (solid), and (iii) almost all 
eigenvalues of the truncated Frobenius-Perron matrix (dashed) . The semiclas- 
sical prediction is mainly influenced from a few eigenvalues of large moduli, 
(b) Comparison of the semiclassical prediction (iii) (dashed) with averaged L 2 
norms of skew Husimi cigcnfunctions (solid) . (c) Logarithmic plot of eigenvalues 
of the truncated Frobenius-Perron matrix (In | A| versus arg A). The peaks in (a) 
are associated to at least one eigenvalue of large modulus. 
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Figure 8: (a) L 2 norms of diagonal Husimi eigenfunctions versus Floquet cigen- 
phases. (b) 530 L 2 norms of most strongly localized skew Husimi eigenfunctions. 
Note the frequent appearance close to the phases of resonances. 
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